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ABSTRACT 

We explore the population of candidate pre-low-mass X-ray binaries in which a neu- 
tron star accretes mass from the wind of a low-mass companion (mass < 2 Mq) in the 
framework of a binary population synthesis study. The simulated accretion-luminosity 
distribution shows a primary peak close to 10^^ erg/s and a secondary peak near 10^^ 
erg/s. The relative contribution of the two peaks depends primarily on the magni- 
tude of the kick velocity imparted to the neutron star at birth. The secondary peak is 
negligible for average kick velocities larger than ^ 200 km/s, but becomes dominant 
for average kick velocities smaller than 50 km/s. Regardless of the relative contri- 
butions of the two peaks, our calculations suggest that pre-low-mass X-ray binaries 
may provide a non-negligible contribution to the population of discrete low-luminosity 
X-ray sources in the Galaxy. 
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1 INTRODUCTION 



Ongoing deep X-ray surveys of selected Galactic regions 
have sparked a renewed interest in low- luminosity Galac- 
tic X-ray sources such as wind-accreting neutron stars. 
The XMM-Newton Galactic Plane Survey (e.g. Warwick 
2002), the XMM-Newton Serendipitous Survey (Watson 
2001; Motch et al. 2002; Motch, Herent & Guillout 2003; 
Watson et al. 2003) and the serendipitous Chandra Multi- 
wavelength Plane Survey (Grindlay et al. 2003) will all probe 
the faint Galactic X-ray point source population. 

The discovery by Chandra of more than 500 previously 
undetected point sources with luminosities above ~ 10^^ 
erg/s in a 2.0° x 0.8° field near the Galactic centre (Wang, 
Gotthelf & Lang 2002) dramatically illustrates the poten- 
tial of such surveys and has already prompted comple- 
menting theoretical studies. Pfahl, Rappaport & Podsiad- 
lowski (2002) considered the population of neutron stars 
that accrete from the wind of an intermediate- or high-mass 
star (mass > 3Mq). They found that these wind-driven 
intermediate-mass and high-mass X-ray binaries may indeed 
account for a dominant fraction of the point source popula- 
tion found by Wang et al. (2002). 

Independently, Bleach (2002) investigated the possibil- 
ity of detecting detached close binaries consisting of a neu- 
tron star and a red dwarf, which may be the progenitors 
of low-mass X-ray binaries (LMXBs). The proposed detec- 
tion mechanism consists of looking for red dwarfs with X-ray 
luminosities larger than those expected from intrinsic X-ray 
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emissions from the star's corona. The excess could then pos- 
sibly be ascribed to the emission of X-rays by a neutron star 
companion accreting mass from the wind of the red dwarf. 

Our aim in this paper is to elaborate on the idea by 
Bleach (2002) and examine the detection probability of pre- 
LMXBs, i.e. detached neutron-star systems with low-mass 
companions (mass < 2Mq), in the context of a binary pop- 
ulation synthesis study. 



2 PRE-LMXBS AND LMXBS 

In close binary systems of stars, mass can be transferred 
from one component to the other when one of the stars loses 
mass in a stellar wind or when one of the stars fills its Roche 
lobe. If the companion star is a compact object which is 
able to accrete some of the transferred mass, the liberation 
of gravitational potential energy gives rise to an accretion 
luminosity 



R 



(1) 



where G is the Newtonian constant of gravitation, Al and R 
are the mass and the radius of the accreting compact object, 
and Mace is the mass-accretion rate. We assume that the 
accretion process yields an X-ray luminosity Lx ~ Lace. 

In the particular case of a wind-accreting neutron star, 
the spin evolution of the neutron star and the possible inhi- 
bition of accretion by strong centrifugal forces may play an 
important role in the study of the accretion process (e.g. 
Pringle & Rees 1972; lUarionov & Sunyaev 1975; Stella, 
White, Rosner 1986; Urpin, Geppert & Konenkov 1998). 
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The treatment of these inhibiting effects is however beyond 
the scope of this initial investigation. Instead, we assume 
that at the onset of mass transfer via the stellar wind, the 
neutron star's radiation pressure and rotation rate is low 
enough for mass accretion to be possible. 

The evolution of binaries containing a neutron star is 
driven by the nuclear evolution of the companion star or by 
the loss of orbital angular momentum via magnetic braking 
and/or gravitational radiation. When the combined effect 
of these processes causes the companion to become larger 
than its Roche-lobe, an X-ray binary may be formed. The 
further evolution of the binary depends critically on the or- 
bital period and the mass of the donor star at the onset of 
Roche-lobe overflow. The formation and evolution of X-ray 
binaries with low-mass donor stars have been discussed ex- 
tensively in the past by, e.g., Verbunt (1993); Verbunt & 
van den Heuvel (1995); Iben & Tutukov (1995); Kalogera & 
Webbink (1996, 1998); Kalogera (1998); Kalogera, Kolb & 
King (1998); and more recently by Podsiadlowski, Rappa- 
port & Pfahl (2002). We therefore refer to these studies and 
references therein for detailed discussions on the various for- 
mation channels and here consider the population of X-ray 
binaries descending from the different channels as a whole. 

For the purpose of this investigation, we define a can- 
didate pre-LMXB as a detached binary in which a neutron 
star accretes mass from the stellar wind of its companion. 
We assume the latter to be on the main sequence, in the 
Hertzsprung gap, or on the giant branch, and to have a 
mass smaller than or equal to 2 Mq . The luminosity asso- 
ciated with the accretion process is required to be larger 
than 10~® Z/0 or ten times the mass donor's coronal X-ray 
luminosity, whichever is larger. We note that our definition 
of pre-LMXBS differs from that of Bleach (2002) in that we 
consider stars with masses up to 2 Mq and binaries that are 
expected to become semi-detached due to the nuclear evo- 
lution of the secondary as well as those that are expected to 
become semi-detached due to angular momentum losses via 
magnetic braking and/or gravitational radiation. 

Similarly, we define a LMXB as a semi-detached system 
in which a neutron star accretes mass from a Roche-lobe 
filling main-sequence, Hertzsprung-gap, or giant-branch star 
with a mass smaller than or equal to 2 Mq . In order to be 
classified as an X-ray binary, we furthermore require the 
mass-accretion rate to be high enough to produce an X-ray 
luminosity in excess of 10~^ Lq. 

For the remainder of the paper, we denote the mass 
and the radius of the neutron star by Mns and -Rns, and 
the mass and the radius of the donor star by Md and Rd, 
respectively. 



3 ASSUMPTIONS AND COMPUTATIONAL 
TECHNIQUE 

We explore the possibility of detecting and identifying 
LMXBs prior to the onset of Roche-lobe overflow from the 
neutron star's companion by means of the BiSEPS binary 
population synthesis code described by WiUems & Kolb 
(2002). The code is based on the analytic approximation 
formulae for the evolution of single stars derived by Hurley, 
Pols & Tout (2000) and follows the main lines of the bi- 
nary evolution algorithm described by Hurley, Tout & Pols 
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Figure 1. The logarithm of the wind mass-loss rates (in 
Af0yr~^) used in our calculations, as a function of the stellar 
mass and radius. The lines labelled ZAMS, TMS, and BGB rep- 
resent the radii of the stars at the zero-age main sequence, the 
terminal main sequence, and the base of the giant branch, respec- 
tively. 



(2002). The binary orbits are assumed to be circular and 
the stellar rotation rates are kept synchronised with the or- 
bital motion at all times. We furthermore restrict ourselves 
to Population I stellar compositions. 

For the purpose of this investigation, the mass-loss rates 
from stellar winds given by Hurley et al. (2000), which are 
limited to evolved stars and massive main-sequence stars, 
must be supplemented by a prescription for the winds from 
low-mass main-sequence stars. Although these winds are not 
fully understood yet, it is generally accepted that they de- 
pend sensitively on both the mass and the age of the donor 
star (e.g. Wood et al. 2002, and references therein). Young 
solar-type stars and M dwarfs, for instance, are thought to 
have mass-loss rates which may be as high as 10~^^ Mq per 
year, while older G dwarfs such as the Sun have mass-loss 
rates of the order of Mq per year (Hartmann 1985, 

Lim & White 1996, van den Oord & Doyle 1997, Wargelin 
& Drake 2001, Sackmann & Boothroyd 2003). In view of 
these still existing uncertainties, we here simply adopt a 
constant mass-loss rate of 10~^^ Mq per year for low-mass 
main-sequence stars (see also Bleach 2002). The resulting 
mass-loss rates are shown in Fig. for donor stars with 
masses up to 2 Mq and evolutionary stages from the zero- 
age main sequence up to the tip of the giant branch. The 
lines labelled ZAMS, TMS, and BGB represent the radii of 
the stars at the zero-age main sequence, the terminal main 
sequence, and the base of the giant branch, respectively. The 
maximum evolutionary age considered for the calculation of 
the mass-loss rates is 10 Gyr, so that only stars with masses 
larger than 1 Mq significantly evolve away from the zero-age 
main sequence. 
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We estimate the mean mass-accretion rate onto a star 
orbiting in the stellar wind of its companion by means of 
a standard Bondi-Hoyle-Lyttleton formalism (Hoyle & Lyt- 
tleton 1941, Bondi & Hoyle 1944). The accretion rate for 
binaries with circular orbits then depends on the orbital 
separation a, the binary mass ratio q, the donor star's ra- 
dius Rd , the wind mass- loss rate Md , and the wind velocity 
Vy,. We parametrise the latter as a fraction /3w of the escape 
velocity Ve = {2GMd/Rd)^^^ at the surface of the mass- 
losing star (see, e.g., Hurley et al. 2002). The mean wind 
mass-accretion rate can then be cast in the form 

In the particular case of a wind-accreting neutron star the 
mass ratio is given by g = Mt^s/Md- 

The detection mechanism for pre-LMXBs proposed by 
Bleach (2002) furthermore requires a prescription for the 
intrinsic X-ray luminosity of the neutron star's companion. 
Since this luminosity is thought to have its origin in dynamo 
processes generated by the interaction between convection 
and rotation, it can be expected to be correlated with the 
star's rotational angular velocity (e.g. Belvedere, Chiuderi 
& Paterno 1982; Noyes et al. 1984). The existence of such a 
correlation was explored extensively from an observational 
point of view by, e.g., Pallavicini et al. (1981); Marilli & 
Catalano (1984); Vilhu (1984); Doyle (1987); Fleming, Gioia 
& Maccacaro (1989); Hempelmann et al. (1995); Randich et 
al. 1996; Singh et al. (1999); Gondoin (1999). The numer- 
ous investigations however show a clear lack of consensus on 
whether or not the correlation exists and, if it exists, what 
the precise relationship is between the rotational angular ve- 
locity and the coronal X-ray luminosity. This deficiency may 
be related to the use of different stellar samples and observa- 
tional selection criteria (Hempelmann et al. 1995, Gondoin 
1999). 

A point of lesser dispute seems to be the existence of 
an upper limit on the intrinsic stellar X-ray luminosity, al- 
though its origin remains to be unambiguously identified 
(Waher 1982, Byrne et al. 1984, Vilhu & Walter 1987, Jar- 
dine &i Unruh 1999, Singh et al. 1999) . The upper limit may 
be related to the finite size of the stellar surface in the sense 
that the magnetic activity can only increase with increasing 
rotational angular velocities until the whole stellar surface 
is covered with magnetically active regions. Based on this 
idea, Fleming et al. (1989) showed that the coronal satura- 
tion luminosity Lcor could be related to the square of the 
stellar radius Rd as 

log Lcor = -2.9 + 2 log i?d, (3) 

where both Lcor and Rd are expressed in solar units (see 
their Fig. 3). 

Rather than face the uncertainties regarding the rela- 
tion between stellar rotation and coronal X-ray activity, we 
here opt to use a prescription for the upper limit on the 
coronal X-ray luminosity. We use the mass of the convec- 
tive envelope, A/conv, to distinguish between coronally active 
and non-active stars and, somewhat arbitrarily, set the limit 
for coronal X-ray activity at Mconv = 0.01 Md. Stars with 
smaller convective envelopes are assumed to have negligible 
coronal X-ray luminosities, while stars with larger convective 
envelopes are assumed to have coronal X-ray luminosities 



given by Eq. (gj. The upper limit may significantly over- 
estimate the coronal X-ray luminosity of giant-type stars 
with slow rotation rates and correspondingly weak dynamos 
(Schroder, Hiinsch & Schmitt 1998; Gondoin 1999). The slow 
rotation, however, only applies to giants in wide binaries 
where tidal forces are negligible. Giants in closer binaries 
are subjected to tidal torques which tend to synchronise the 
giant's rotation with the orbital motion of the companion. 
A significant dynamo action may therefore be sustained in 
these stars during much more evolved stages of stellar evolu- 
tion than for single stars or stars in wide binaries. Dempsey 
et al. (1993, 1997) have shown that in the particular cases of 
RS Canum Venaticorum (RS CVn) and BY Draconis (BY 
Drac) binaries the coronal X-ray luminosity of the giant- 
type component may still reach values up to 0.1 Lq. To be 
on the safe side, we therefore still use the upper limit given 
by Eq. @ for stars that have evolved beyond the main- 
sequence. Where relevant, we will discuss the possible im- 
plications of this on our results. In any case, we expect the 
consistent use of the upper limit for the coronal X-ray lumi- 
nosity to imply that our selection criteria for the identifica- 
tion of pre-LMXBs are systematically too strict. 



4 POPULATION SYNTHESIS 

We evolved a large number of binaries with initial compo- 
nent masses ranging from 0.1 to 60 Mq and with initial or- 
bital periods between 0.1 to 10 000 days. We used logarithmi- 
cally spaced grids consisting of 60 grid points for the initial 
stellar masses and 200 grid points for the initial orbital peri- 
ods. The maximum evolutionary age considered was 10 Gyr. 
If a binary undergoes an asymmetric supernova explosion 
leading to the birth of a neutron star, 200 random kick ve- 
locities are generated from a Maxwellian velocity distribu- 
tion for the magnitude of the kick velocities. The evolution 
of the binaries surviving the explosion is then followed until 
the imposed age limit of 10 Gyr. 

The outcome of the binary evolution calculations de- 
pends on the parameters adopted for the various evolution- 
ary processes a binary may be subjected to. We therefore 
vary the input parameters one at a time and compare the 
resulting populations of pre-LMXBs with that of a refer- 
ence model. For the latter, hereafter referred to as model A, 
we adopt the same set of input parameters as in Willems 
& Kolb (2002), i.e. the common-envelope ejection efficiency 
parameter takes the value acE = 1.0, and the Maxwellian 
kick- velocity dispersion takes the value Ukick = 190 km/s. In 
addition, we here adopt a standard value of 1.0 for the wind- 
velocity parameter /3w The different parameters adopted in 
the other population synthesis models are summarised in Ta- 
ble0 In the case of the different supernova-kick models, the 
average kick velocities i?kick adopted in the Maxwellian ve- 
locity distributions are listed in addition to the velocity dis- 
persions akick- We also note that the variations in the wind- 
velocity parameter /3w considered in models MAI and MA2 
are equivalent to considering variations in the wind mass-loss 
rate Md, but that the parameter /3w is a more appropriate 
population synthesis parameter since it has a greater impact 
on the determination of the mean mass-accretion rate. This 
is readily illustrated by means of Eq. ^ which shows that, 
to a first approximation, an increase in /3w by a factor of 2 
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corresponds to a decrease in Md by a factor of the order of 
10. Models CORl, COR2, and COR3 are introduced to as- 
sess the sensitivity of our results to the prescription for the 
coronal X-ray luminosity. In particular, the use of a constant 
value for L^or in models CORl and COR2 allows us to ac- 
count for the possible effects of overestimating the coronal 
X-ray luminosity of giant-type stars by the use of Eq. @. 

The contribution of a system to the population of pre- 
LMXBs depends on the probability distribution functions 
of its initial parameters and on the time the system spends 
as a member of the population. We assume the initial mass 
Afi of the primary to be distributed according to the initial 
mass function 

Mi/Mq < 0.1, 

C(Mi) = { 0.384:15 M^^ 0.1 < M1/M0 < 0.75, (4) 
0.23556 M"^ '^ 0.75 < Mi/Mq < 00, 

the initial mass ratio q = M2/M1 according to 



n(q) = 



1 < g < 1, 
,0 q>l, 

and the initial orbital separation a according to 







a/RQ < 3 or a/i?© > lO", 



0.12328 fl-i 3 < q/Rq < 10* 



(5) 



(6) 



(see Willems & Kolb 2002, and references therein). 

Each time a binary from our initial grid of parameters 
evolves into a pre-LMXB, its weighted contribution is added 
to the probability density function (PDF) describing the dis- 
tribution of newborn pre-LMXBs in the (t, Md, Porb, ^acc)- 
space at the time t = tb of its formation. The finite lifetime 
of the pre-LMXBs is taken into account by subtracting the 
same weighted contribution from the distribution function 
at the time t — td when the system ceases to be a pre-LMXB. 
In doing so, we implicitly assume that the stellar and orbital 
parameters of a pre-LMXB do not change significantly once 
it is formed. We will discuss the circumstances under which 
this assumption breaks down and the possible implications 
thereof in the final section of the paper. For now, we remark 
that this assumption does not affect the determinations of 
the birthrates and the total number of systems obtained af- 
ter integration over all system parameters. 

The resulting distribution function for the population of 
pre-LMXBs is subsequently convolved with a constant star- 
formation rate and normalised so that the integral over all 
systems found is equal to one. 



5 DONOR STAR MASSES AND ORBITAL 
PERIODS 

The normalised distribution function of pre-LMXBs in the 
(log Md, log Porb)-plane is displayed in Fig. |5|for the most 
illustrative of the population synthesis models listed in Ta- 
ble Bins containing systems that actually evolve into 
LMXBs within the imposed age limit of 10 Gyr are outlined 
with a black solid line. 

There are three reasons why systems satisfying our pre- 
LMXB criteria do not show up in the corresponding LMXB 



population. Donor stars with mass Md ^ 1 Mq evolve too 
slowly to fill their Roche lobe as a result of their nuclear 
evolution within the imposed age limit of 10 Gyr. They 
can therefore only become LMXBs if their orbital period is 
shorter than the bifurcation period of ~ 1 day, which sepa- 
rates diverging from converging systems (Pylyser & Savonije 
1988, 1989). Donor stars with mass Md > 1.4 Mq in systems 
with Porb ^2.5 days on the other hand may fill their Roche 
lobe when the donor star ascends the first giant branch, but 
the resulting mass transfer is dynamically unstable so that 
a common-envelope phase ensues instead of a LMXB phase. 
Systems with Porb ^ 800 days, finally, are too wide for the 
donor star to fill its Roche lobe prior to the AGB stage. 

The lack of pre-LMXBs in the upper-left corner of each 
of the plots presented in Fig. |5| is related to our treatment 
of the coronal X-ray luminosity of the donor stars. Since 
stars less massive than ~ 1 Mq are essentially unevolved 
and have appreciable convective envelopes, the adopted up- 
per limit for the coronal X-ray luminosity is proportional 
to the square of the zero-age main sequence radii [Eq. Q]. 
The decrease of the latter with decreasing mass then implies 
that binaries with lower-mass donor stars can satisfy the cri- 
terion Lace > 10 I/cor Up to longer orbital periods. The exact 
position of the lower border of the wedge in the upper-left 
corner of the figures introduced by this criterion depends 
on the wind-velocity parameter /3w and on the assumptions 
made for the determination of the coronal X-ray luminosity 
Lcor. Increasing the wind- velocity parameter (model MA2) 
yields lower mass- accretion rates and thus lower accretion lu- 
minosities. Correspondingly, the bottom of the wedge shifts 
to lower orbital periods. Adopting a constant upper limit for 
the coronal X-ray luminosity (models CORl and COR2, not 
shown) instead of Eq. ^ renders the border of the wedge 
introduced by the Lace > 10 Lcor criterion independent of 
the mass of the donor star. 

Stars more massive than ~ IMq, on the other hand, 
may evolve away from the ZAMS within the imposed age 
limit of 10 Gyr. However, the deep convective envelopes and 
large radii developed on the giant branch yield high up- 
per limits for the coronal X-ray luminosity so that only 
few wind-accreting neutron stars manage to top this limit 
with their accretion luminosity. Most of the pre-LMXBs with 
Md > 1 Mq therefore also have main-sequence donor stars. 
The rather sharp vertical border of the wedge in the upper- 
left corner of the figures at Md ~ 1 Mq is then related to 
the assumed lower limit Afconv = 0.01 Add for the mass of the 
convective envelope of coronally active stars. Main-sequence 
stars more massive than ~ 1 Mq have smaller or no con- 
vective envelopes and are thus assumed to have negligible 
coronal X-ray luminosities. The pre-LMXB luminosity cri- 
terion therefore reduces to Lace > 10~® Lq, which is easily 
satisfied. If the minimum convective envelope mass neces- 
sary for X-ray activity is increased to, e.g., Meonv = 0.1 Afd 
(model C0R3, not shown), the sharp vertical border of the 
wedge shifts to Md ~ 0.8 Mq. 

The group of systems with orbital periods longer than 
~ 1000 days and donor star masses between 1.0 and 1.5 M© 
in the low kick-velocity dispersion models shown in Fig. |5| 
corresponds to systems with initial orbital separations wide 
enough to avoid any kind of Roche-lobe overflow prior to 
the formation of the neutron star. These systems easily sur- 
vive the supernova explosion of the neutron star's progenitor 
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Table 1. Population synthesis model parameters. 



model 




/3w 


""kick 


fkick 


-^cor/ -^'O 


Mconv/Md 


A 


1.0 


1.0 


190 km/s 


330 km/s 


Eq. 


0.01 


KO 


1.0 


1.0 


no kicks 




Eq. 


0.01 


KM25 


1.0 


1.0 


25 km/s 


43 km/s 


Eq. 13 


0.01 


KM50 


1.0 


1.0 


50 km/s 


87 km/s 


Eq. 13 


0.01 


KM75 


1.0 


1.0 


75 km/s 


130 km/s 


Eq. 13 


0.01 


KMIOO 


1.0 


1.0 


100 km/s 


173 km/s 


Eq. 13 


0.01 


KM300 


1.0 


1.0 


300 km/s 


520 km/s 


Eq. 13 


0.01 


KM400 


1.0 


1.0 


400 km/s 


693 km/s 


Eq. 13 


0.01 


CEl 


0.2 


1.0 


190 km/s 


330 km/s 


Eq. 13 


0.01 


CE2 


5.0 


1.0 


190 km/s 


330 km/s 


Eq. 13 


0.01 


MAI 


1.0 


0.5 


190 km/s 


330 km/s 


Eq. 13 


0.01 


MA2 


1.0 


2.0 


190 km/s 


330 km/s 


Eq. 13 


0.01 


CORl 


1.0 


1.0 


190 km/s 


330 km/s 


10"'^ (constant) 


0.01 


COR2 


1.0 


1.0 


190 km/s 


330 km/s 


10~* (constant) 


0.01 


COR3 


1.0 


1.0 


190 km/s 


330 km/s 


Eq. 13 


0.1 



when the average kick velocity is small compared to the rel- 
ative orbital velocity of the component stars, but generally 
get disrupted for higher average kick velocities. 

The narrow horizontal ridge of systems near Pmh ~ 10 
days in the low kick-velocity dispersion models ((Tkick ^ 100 
km/s) is associated with systems that spiralled-in during 
a common-envelope phase triggered by Roche-lobe overflow 
from the neutron star's progenitor prior to its supernova ex- 
plosion. The range of orbital periods available for this evo- 
lutionary channel is limited by the requirements that the 
pre-common-envelope orbit must be small enough to initi- 
ate Roche-lobe overflow prior to the supernova explosion and 
wide enough to provide enough orbital energy to expel the 
common envelope before the core of the Roche-lobe over- 
flowing star merges with the spiralling-in companion. After 
the expulsion of the envelope, the change in the orbital con- 
flguration caused by the supernova explosion depends on the 
amount of mass lost from the system and on the kick veloc- 
ity imparted to the neutron star at birth. For a detailed 
account of the effects of asymmetric supernova explosions 
on the orbital characteristics of binary stars, we refer to 
Kalogera (1996). 

In the case of small supernova kick- velocity dispersions, 
the effect of the mass loss dominates over the effect of the 
kick. The small spread of orbital periods available for the 
common-envelope scenario is then preserved by the one-to- 
one relation between the post-supernova orbital period and 
the mass lost from the system (e.g. Verbunt 1993). Further- 
more, since binaries in which more than half of the total 
mass is lost from the system become unbound after a sym- 
metric supernova explosion, the probability of survival is 
larger for binaries with more massive donor stars. 

In the case of higher supernova kick- velocity dispersions 
the effect of the kick dominates over the effect of the mass 
loss. This facilitates the survival of binaries with lower-mass 
donor stars and increases the spread of the post-supernova 
orbital periods. The horizontal ridge around Porb ~ 10 days 
therefore widens and shifts to lower donor star masses. 

The common-envelope ejection efflciency parameter 
QCE only affects systems with orbital periods less than 
~ 1000 days. A decrease of acE (model CEl) yields smaller 
orbital separations at the end of the common-envelope phase 



and thus shifts the horizontal ridge near Porb ~ 10 days to 
lower orbital periods. This effectively widens the gap be- 
tween the long- and short-period systems. Increasing the 
common-envelope efficiency acE has the opposite effect. 

We conclude this section by noting that the position 
of the high-density areas in Fig. |21 also depends on the 
adopted initial mass-ratio distribution. Mass ratio distribu- 
tions favouring lower initial mass ratios (e.g. n{q) — 1/q) 
shift the high-density areas to lower donor star masses, while 
mass ratio distributions favouring higher initial mass ratios 
(e.g. n{q) = q) shift the high-density areas to higher donor 
star masses. 



6 X-RAY LUMINOSITIES 

The normalised distribution functions of pre-LMXBs in the 
(log Porb, log Lace)- and the (log Md, log Lace )-plane are dis- 
played in Figs. 13 and 21 respectively. Since the qualitative 
differences between the models listed in Tableware small, 
we limit the figures to models KM50, A, and MA2. 

The bulk of the systems shown in Fig. |3 have loga- 
rithmic accretion luminosities that are linearly correlated 
with the logarithm of the orbital period. The slope of -4/3 
arises from Kepler's third law and the factor l/a^ in Eq. 
12)l for the mean mass-accretion rate Mace- Since most of 
the systems following the correlation have main-sequence 
donor stars, the small spread around the correlation is the 
result of the limited range of donor star masses and radii 
and the adopted constant wind mass-loss rate of 10~^^ Mq 
per year [see Eq. In the case of low supernova kick- 

velocity dispersions, the two distinct higher-density regions 
along the Lace oc P^^^^^ relation correspond to the two 
groups of systems identified in the low kick- velocity disper- 
sion models shown in Fig.|21 The high accretion luminosities 
around 10^^ erg/s stem from systems on the horizontal ridge 
around Porb ~ 10 days, while the lower accretion luminosi- 
ties around 10^* erg/s stem from the long-period systems 
with Porb ^ 1000 days. For higher kick-velocity dispersions, 
only the first group of systems survives. The accretion lumi- 
nosities of the pre-LMXBs on the Lace oc Po^b''^ furthermore 
decrease with increasing values of the wind-velocity param- 
eter /3w 
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Figure 2. Normalised distribution of pre-LMXBs in the (log Mj, log Porb)-pla'ne. The outlined boxes correspond to bins containing 
pre-LMXBs which actually evolve into LMXBs within the imposed age limit of 10 Gyr. 



The branch of systems located to the right of the 
Lace oc P~J^^^ correlation corresponds to pre-LMXBs with 
giant branch donor stars. Most of these already appeared 
as pre-LMXBs when the donor star was still on the main- 
sequence, but the growth of the convective envelope and the 
associated increase in the coronal X-ray activity during the 
post-main-sequence evolution temporarily concealed the X- 
ray accretion luminosity emitted by the wind-accreting neu- 
tron star. The functional form of the mean mass-accretion 
rate Mace [Eq. lO] implies that systems with longer orbital 



periods require more evolved donor stars, i.e. stars with 
larger radii and higher wind mass-loss rates, to satisfy the 
imposed Lace > lOLcor criterion. However, the larger radii 
in turn imply a larger upper limit for the coronal X-ray 
luminosity Leor [Eq. so that at the time of reappear- 
ance the 10 Leor threshold itself will be higher for systems 
with wider orbital separations. Hence, the accretion luminos- 
ity at which the systems reappear as pre-LMXBs increases 
with increasing orbital periods. Their reappearance also de- 
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Figure 3. Normalised distribution of pre-LMXBs in the (log Porb, log Lacc)-plane. The outlined boxes correspond to bins containing 
pre-LMXBs which actually evolve into LMXBs within the imposed age limit of lOGyr. 



pends critically on whether or not the accretion luminosity 
increases faster than the coronal X-ray luminosity. 

The position of the giant-branch stars in the 
(log Porb, log Lacc)-plane clearly depends on the assumptions 
adopted for the determination of Lace and L^or ■ Since higher 
wind velocities yield lower mass-accretion rates, an increase 
in the wind- velocity parameter /3w as in model MA2 requires 
even more evolved donor stars to satisfy the Lace > lOLcor 
criterion than in model A. The reappearance threshold of 
these systems consequently moves to higher accretion lumi- 
nosities. In the case of a constant coronal X-ray luminosity 
(models CORl and COR2, not shown), the accretion lumi- 
nosity at which systems with giant branch donor stars reap- 
pear as pre-LMXBs is independent of the orbital period. The 
temporary concealment of the accretion luminosity proba- 
bly does not occur for the longer period systems (Porb ^ 100 
days) in which tidal forces are too weak to keep the donor 
star's rotation synchronised with the orbital motion of the 
companion. The decrease of the donor star's rotational angu- 
lar velocity during its post-main-sequence evolution is here 
more likely to be accompanied by a decrease in the coronal 
X-ray activity so that the wind-accretion luminosity in these 
binaries has no problem topping the coronal X-ray luminos- 
ity at any time after they first appear as a pre-LMXB. 

The linear correlation between log Lace and log Porb for 
the bulk of the systems presented in Fig.Oessentially makes 
Fig.|l]a mirror image of Fig.|5| The low kick- velocity disper- 
sion models hence again show two distinct groups of systems 
at high and low accretion luminosities, which are associated 
with the horizontal ridge and the long-period systems ob- 
served in Fig. 121 In our standard model (model A), even 
systems with donor stars less massive than 1.0 are seen 
to give rise to accretion luminosities up to lO'^'^ erg/s. An 
increase in the wind- velocity parameter /3w shifts all systems 
to lower accretion luminosities. 

Finally, the one-dimensional distribution functions of 
the pre-LMXB accretion luminosities are presented in Fig.|3 
for the same population synthesis models as in Fig. |21 
We distinguish between systems that actually evolve into 
LMXBs within the imposed age limit of 10 Gyr and those 



that do not by means of Ught and dark grey shadings, respec- 
tively. The distribution of pre-LMXB accretion luminosities 
in our standard population synthesis model shows a primary 
peak around 10^^ erg/s and a small secondary peak around 
10^* erg/s. The latter becomes more important with decreas- 
ing kick-velocity dispersions and dominates when (Tkick ~ 25 
km/s. This shift towards lower accretion luminosities is re- 
lated to the increasing contribution of the long-period sys- 
tems with Porb > 1000 days. The width of the peak at 10^^ 
erg/s depends primarily on the common-envelope ejection 
efficiency parameter qce. The peak becomes narrower with 
decreasing values of acE due to the shrinking range of post- 
supernova orbital periods around Porb ~ 10 days. Increas- 
ing the wind- velocity parameter /3w moves the distribution 
function to lower X-ray luminosities, while a decrease in /3w 
has the opposite effect. A similar tendency was found by 
Pfahl et al. (2002). It is also interesting to note that the 
peak luminosity of 10"^^ erg/s is comparable to the lumi- 
nosity originating from non-magnetic white dwarfs accret- 
ing mass from a Roche-lobe filling companion (cataclysmic 
variables). Such systems may be distinguishable from pre- 
LMXBs based on the harder X-ray spectrum expected for 
wind-accreting neutron stars. In addition, white dwarfs can 
possibly be detected in the optical or UV wavelengths, pro- 
vided that they are bright enough. 



7 RELATIVE AND ABSOLUTE NUMBER OF 
PRE-LMXBS IN THE GALAXY 

The relative number of pre-LMXBs evolving into LMXBs 
within the imposed age limit of 10 Gyr are listed in the first 
column of Table |5| In general, about 20 to 30 per cent of 
all systems satisfying our pre-LMXB criteria evolve into a 
LMXB in a time span of less than 10 Gyr. This fraction 
decreases significantly for kick-velocity dispersions smaller 
than ~ 50 km/s (model KM50) due to the growing contri- 
bution of systems with orbital periods longer than ~ 1000 
days. A decrease of the common-envelope efficiency parame- 
ter acE (model CEl) on the other hand increases the relative 
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Figure 4. Normalised distribution of pre-LMXBs in tlie (log Md, log Lace )-plane. The outlined boxes correspond to bins containing 
pre-LMXBs which actually evolve into LMXBs within the imposed age limit of lOGyr. 



Table 2. Relative number of pre-LMXBs evolving into LMXBs 
within the imposed age limit of 10 Gyr and relative number of 
LMXBs descending from a pre-LMXB. 



pre-LMXBs evolving LMXBs descending 



model 


into LMXBs 


from pre-LMXBs 


A 


31.9% 


43.5% 


KO 


21.3% 


0.1% 


KM25 


12.7% 


14.3% 


KM50 


16.1% 


10.8% 


KM75 


24.1% 


8.5% 


KMIOO 


28.7% 


9.7% 


KM300 


30.1% 


48.7% 


KM400 


29.5% 


48.2% 


CEl 


40.7% 


14.9% 


CE2 


33.2% 


73.6% 


MAI 


26.1% 


36.9% 


MA2 


41.2% 


39.1% 


CORl 


45.8% 


43.2% 


COR2 


24.7% 


43.5% 


COR3 


26.5% 


43.5% 



range 30 to 50 per cent. The LMXBs that did not evolve 
through a pre-LMXB phase all descend from intermediate 
mass X-ray binaries (IMXB) with initial donor star masses 
outside the mass range considered in this investigation (see 
also King & Ritter 1999; Podsiadlowski & Rappaport 2000; 
Kolb et al. 2000; Pfahl, Rappaport & Podsiadlowski 2003). 
A substantial fraction of LMXBs thus evolves through an X- 
ray bright wind-accretion phase prior to the onset of Roche- 
lobe overflow from the neutron star's companion. 

For conclusion, we estimate the order of magnitude of 
the birthrate and the absolute number of pre-LMXBs cur- 
rently populating the Galaxy. To this end, we assume all 
stars to be in binaries and we normalise the star formation 
rate to produce a Galactic type II supernova rate of 10~^ 
per year (Cappellaro, Evans & Turatto 1999). The result- 
ing birthrates and absolute numbers are listed in Table |21 
Based on these calculations, we expect the current popula- 
tion of pre-LMXBs in the Galactic disk to exist of about 10*- 
10^ systems, so that they may account for a non-negligible 
part of the X-ray output in large scale surveys performed by 
Chandra and XMM. 



number of pre-LMXBs evolving into LMXBs since it boosts 
the contribution of the short-period systems. Similarly, in- 
creasing the wind-velocity parameter /3w (model MA2) de- 
creases the accretion luminosities of the longer period sys- 
tems below the Lace > 10 Leor threshold so that again the 
relative number of pre-LMXBs evolving into LMXBs in- 
creases. If Eq. is replaced by a constant coronal X-ray 
luminosity of 10~^ Lq (model CORl), a large fraction of the 
low-mass systems with Md < 1 Mq can no longer satisfy the 
iace > 10 Leor Criterion. The disappearance of these systems 
causes an increase of the relative number of pre-LMXBs that 
will evolve into LMXBs. 

Although a detailed population synthesis study of the 
formation and evolution of LMXBs is beyond the scope if 
this investigation, it is interesting to consider the relative 
number of LMXBs that are expected to evolve through a 
phase satisfying our pre-LMXB criteria. The relative num- 
bers are listed in the second column of Table|5|and typically 



8 DISCUSSION 

We have investigated the population of candidate pre-low- 
mass X-ray binaries according to the detection mechanism 
proposed by Bleach (2002). The mechanism is based on the 
search for low-mass stars with X-ray luminosities that are 
too large to be attributed to the stellar coronal activity, 
and on the subsequent identiflcation of neutron star binaries 
amongst these stars. The X-ray excess in these binaries can 
then possibly be ascribed to the accretion of mass by the 
neutron star from the stellar wind of its companion. 

In our study, we used the BiSEPS binary popula- 
tion synthesis code described by Willems & Kolb (2002). 
We limited ourselves to neutron star binaries with main- 
sequence, Hertzsprung-gap, or giant-branch donor stars of 
mass Md < 2 Mq and simulated the X-ray excess by the re- 
quirement that the X-ray luminosity generated by the wind- 
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Figure 5. Normalised distribution functions of pre-LMXB accretion luminosities. Light and dark grey shadings correspond to the 
contributions of systems that do and do not evolve into LMXBs within the imposed age limit of 10 Gyr, respectively. 



accreting neutron star be at least 10 times larger than the 
coronal X-ray luminosity of the mass-losing star. The X-ray 
luminosity associated with the accretion process was deter- 
mined under the assumptions that no accretion-inhibiting 
effects such as strong centrifugal forces are present and that 
all gravitational potential energy of the infalling matter is 
converted into X-ray radiation. The maxiumum coronal X- 
ray luminosity for a star of a given mass and radius was 
assumed to be proportional to the square of the star's ra- 
dius. 

It follows that, in the case of low supernova kick- 
velocity dispersions, the population of pre-LMXBs can be 



divided into two groups of systems: a long-period group with 
forb ^ 1000 days and a short-period group with Porb ~ 10 
days. The division yields a bimodal distribution for the ac- 
cretion luminosities of the candidate pre-LMXBs. The long- 
period systems produce a peak in the distribution around 
lO'^* crg/s, while the short-period systems produce a peak 
near 10"^^ erg/s. For kick- velocity disperions larger than 
~ 100 km/s only the peak near 10"'^ erg/s remains. The 
luminosity at which the distribution peats furthermore de- 
pends sensitively on the velocity of the donor star's wind. 

These results are derived under the assumption that the 
stellar and orbital parameters, and in particular the orbital 
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Table 3. Birthrates and total number of Galactic pre-LMXBs 
for each of the population synthesis models listed in Table 111 The 
models are normalised to a Galactic type II supernova rate of 
10~^ per year. 



model birthrate (yr ^) total number 
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period and the accretion luminosity, of pre-LMXBs remain 
unchanged after their formation. This assumption breaks 
down for binaries with Afd ^ 1 Mq and Porb ^ 1 day in which 
the evolution is governed by the loss of angular momentum 
via gravitational radiation and/or magnetic braking, and for 
binaries with Md ^ 1 Mq in which the evolution is governed 
by the nuclear evolution of the donor star. 

First, the decrease of the orbital period resulting from 
the angular momentum loss in binaries with Md < 1 Mq and 
Porb ^ 1 day yields an increase of the mean mass-accretion 
rate and thus of the accretion luminosity generated by the 
neutron star. However, from Eqs. Q and ||5J it follows that 
an increase of the accretion luminosity by a factor of 10 re- 
quires a decrease of the orbital period by a factor of ~ 5.6. 
For a pre-LMXB with an orbital period of ~ 1 day, this 
comes down to a decrease of the period to ~ 4.3 hours. 
Since neutron star binaries with secondaries more massive 
than ~ 0.5 Mq becomes semi-detached before they can reach 
such short periods, we can expect that the neglect of the or- 
bital evolution does not affect our estimates for the accretion 
luminosities by more than an order of magnitude. 

Second, the nuclear evolution of donor stars more mas- 
sive than ~ 1 Mq may yield a significant increase of the 
star's radius and wind mass-loss rate within the age of the 
Galaxy. The effect of both these evolutionary aspects is 
to increase the mean mass-accretion rate of the neutron 
star. However, for main-sequence stars less massive than 
2 Mq , the increase of the stellar radius increases the mean 
mass- accretion rate by less than an order of magnitude [see 
Eqs. Q and The overall effect on our results can there- 
fore again be expected to be small. The main uncertainty for 
main-sequence donor stars consequently stems from their 
poorly understood wind mass-loss rates (e.g. Wargelin & 
Drake 2001). Pre-LMXBs with donor stars that have evolved 
beyond the main-sequence, on the other hand, have better 
understood mass-loss rates, but contribute too little to the 
total population of newborn pre-LMXBs to have a signifi- 
cant effect on the distribution functions presented. In addi- 
tion, systems with an orbital period of the order of a few 
days will evolve into LMXBs shortly after the secondary 



leaves the main sequence. A significant increase in the X- 
ray accretion luminosity is therefore only possible for wider 
systems which do not initiate Roche-lobe overfiow until the 
secondary is well on its way towards the tip of the giant 
branch or for systems which do not initiate Roche-lobe over- 
flow at all. For these systems, we may underestimate the 
maximum accretion luminosity by several orders of magni- 
tude, but since the life time of a star on the giant branch is 
much shorter than its main sequence life time, the contribu- 
tion of the systems to the pre-LMXB population during this 
phase is smaller than their contribution when the donor star 
was still on the main sequence. In addition, the wider orbital 
separations required for this scenario somewhat counteract 
the increase in the accretion luminosity associated with the 
evolution of the donor star. 

Finally, we estimated the total number of pre-LMXBs 
currently populating the Galaxy to be of the order of lO"- 
10^ For comparison, the estimated number of RS CVn type 
binaries in the Galaxy is of the order of lO" (Watson 1990) 
and the theoretically predicted number of cataclysmic vari- 
ables in the Galaxy is of the order of 10® — 10^ (Kolb & 
Willems 2003, in preparation). Both these types of systems 
have X-ray luminosities comparable to the 10^^ erg/s peak 
luminosity found in our models for candidate pre-LMXBs. 
In addition, the distribution of X-ray luminosities at 10^* 
erg/s may be expected to be dominated by coronally ac- 
tive stars. Our predicted number of pre-LMXBs satisfying 
the X-ray excess criterion is therefore small in comparison to 
the number of these other low-luminosity X-ray sources, but 
nevertheless high enough to possibly account for an inter- 
esting fraction of the low-luminosity discrete X-ray sources 
observed by Chandra and XMM. A systematic search for 
candidate pre-LMXBs based on the X-ray excess criterion 
can in principle be initiated by comparing exisiting X-ray 
(ROSAT, Einstein, Chandra and XMM) and optical (Hip- 
parcos, Tycho, HST GSC) source catalogs. Simultaneous X- 
ray and optical observations by XMM provide an attractive 
alternative. 
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